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Dynamics of the polymorphic A<->B transitions in DNA is compared for two polypurine se- 
quences, poly (dA). poly (dT) and poly (dG). poly (dC), long known to exhibit contrasting properties 
in experiments. In free molecular dynamics simulations reversible transitions are induced by chang- 
ing the size of a water drop around DNA neutralized by Na^ ions. In poly (dG). poly (dC) the 
B<->A transitions are easy, smooth and perfectly reversible. In contrast, a B->A transition in 
poly (dA). poly (dT) dodecamer fragment could not be obtained even though its A-form is stable 
under low hydration. Normal B->A transitions are observed, however, in long poly(dA).poly(dT) 
stretches flanked by GC pairs. An intermediate range of hydration numbers is identifled where oppo- 
site transitions are observed in the two dodecamer fragments, namely, A->B in poly (dA). poly (dT) 
and B->A in poly(dG).poly(dC). With hydration numbers close to the stability limit of the B-form, 
the two sequences exhibit qualitatively different counterion distributions, with a characteristic accu- 
mulation of Na"*" ions next to the opening of the minor groove in poly (dA). poly (dT). This difference 
can explain an increased persistence of poly (dA). poly (dT) DNA towards A-form in crystalline and 
amorphous fibers as compared to solution conditions. The good overall agreement with experi- 
mental data corroborates the general role of the electrostatic condensation mechanism in the A/B 
polymorphism in DNA. 



Introduction 

Reversible transitions between the A and B forms of 
DNAi represent one of the modes for governing protein- 
DNA interactions, with the B-form as the dominant 
biological conformation and the A-form as a high en- 
ergy state adopted temporarily for different purposesSi^. 
These transitions can be induced in several ways in 
vitroi'^'^ and they were extensively studied in the re- 
cent decades. However, a number of interesting is- 
sues are still unclear including the underlying molecu- 
lar mechanisms and the driving forces. Different the- 
oretical approaches were applied in this field during 
its long historjiSiSiifliiLiS. In the recent years realistic 
molecular dynamics (MD) simulations of A^B transi- 
tions became possible including relevant expHcit solvent 
environmen^iiii*iS*iS*ii'i£*ifi*21i. Despite the limited du- 
ration of trajectories and some force field artifactsiSl, 
these studies helped to get new insights in the putative 
molecular mechanisms of the DNA A/B polymorphism. 
A new approach recently found in this field consists in 
using small water drops as solvent shells around DNA 
fragments22,. By changing the size of the drop, reversible 
A^B transitions can be induced, with hydration num- 
bers in rough agreement with experimental values-". 

The water drop simulation approach is the first and 
the only presently known technique to obtain reversible 
B^A transitions in DNA in silico. The established 
mechanism of such transitions essentially consist in an 
intra-duplex electrostatic condensation of the two DNA 
strands. With the drop size gradually reduced, the equi- 
librium is shifted towards the A-form due to inversion of 
electrostatic interactions across the major DNA groove 
when the local cation concentration exceeds a certain 



transition level. The important role of this particu- 
lar type of interactions is well established in the ear- 
lier literature^SiSi. Similar B^A transitions were ob- 
served experimentally in poly(dG).poly(dC) fragments 
with multivalent cations bound in the major groove^^, 
which could be reproduced in simulationaifii. This mech- 
anism is perhaps generally responsible for the A/B poly- 
morphism, but a number of other ideas concerning its 
nature are also being discussed. Among them one can 
mention the factor of economy of hydration of the sugar 
phosphate backboneS^, the base stacking interactions^ 
and the hydrophobic effecfeiii^S. It is interesting, there- 
fore, to know if the A^B transitions in water drop simu- 
lations that are relatively simple to study and rationalize 
are really similar to those in experimental conditions. To 
this end we should check if different aspects of the A/B 
polymorphism in silico correspond to the experimental 
results for the most part obtained for DNA in fiber crys- 
tals under low relative humidity and in concentrated so- 
lutions of some non-polar solvents. Among such aspects 
the sequence dependent nature of the A/B polymorphism 
is perhaps most intriguing and well documented. 

Sequence effects upon the B— >A transition attract 
long-standing interest. In fiber crystals, DNAs with 
different regular sequences often exhibit dissimilar ca- 
pacities to adopt the two forms^i^S. The same trends 
are observed in solution, moreover, here the pattern re- 
vealed is much richer because a large variety of sequences 
could be probed, with their propensities towards A- and 
B-forms measured quantitativelji^. It was shown that, 
during transitions in long natural DNA, different re- 
gions switch between the two forms in a complex se- 
quence specific ordes^i. The A/B-philicity of different 
base pair steps and triplets was examined systematically 
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in synthetic fragments. It could be parametrized and 
used successfully for predicting the properties of mixed 
sequence^SiSyi. All these effects certainly originate from 
interactions between stacked bases and some of them 
could be rationalizedSiiSiii. 

The measured A/B propensities appear such that 
consecutive base pair steps tend to compensate one 
another, therefore, generic DNA sequences commonly 
are neither A- nor B-philio21. A strong differ- 
ence is observed, however, between the stretches of 
consecutive guanines and adenines (G-tracts and A- 
tracts, respectively). Double helical DNA fragments 
with these two sequences exhibit contrasting physico- 
chemical properties as established by a number of differ- 
ent m:th:dj"''' '^ -'' '^ -''"^ -'^ ''^ '" ^" ^^ A 
systematic review of these data has been published 
recently^. In G-tracts, the B— >A transition is very 
easy. The poly(dG).poly(dC) DNA is generally prone 
to adopt the A-DNA conformation in conditions where 
random sequences stay firmly in the B-formSi2£iifi. In 
contrast, in A-tracts, the B^A transition is particularly 
difficult^ Notably, a B^A transition was never ob- 
served in fiber crystaUine poly (dA). poly (dT) because in 
standard transition conditions it is reportedly prone to 
denaturation or alternative transformation a^i^^ . Solvent 
cations were found to bind in the major grooves of G- 
tracts and in the minor grooves of A-tracts and not vice 
versa, which is probably related and perhaps responsible 
for the A/B propensities of these sequences^S*^. This re- 
markable physico-chemical difference is particularly im- 
portant in the context of the intriguing roles of A-tracts 
and G-tracts in the genome functioning. Various A-tract 
sequences are overrepresented in both procariotic and eu- 
cariotic genomes^. They are long known to cause DNA 
curvature in solution that can affect binding to proteins 
and formation of nucleosomes^S*^. Long G-tracts are 
also widespread in eucariotic as well as archaebacteria 
genomes^. 

Here we report about the results of systematic sim- 
ulation studies of A<-^B transitions in a series of DNA 
fragments with these two characteristic sequence motifs. 
The established experimental preferences of A- and G- 
tracts towards B and A forms, respectively, are quali- 
tatively well reproduced. In poly(dG).poly(dC) dode- 
camer, the B^A transition is easy, smooth and per- 
fectly reversible. In contrast, a B^A transition in similar 
poly (d A) .poly (dT) fragment could not be obtained even 
though its A-form was found stable under low hydration. 
At the same time, normal B— s-A transitions are observed 
in long poly (dA). poly (dT) stretches flanked by GC pairs. 
An intermediate range of hydration numbers is found 
where opposite transitions are observed in the two do- 
decamer fragments, namely, A^B in poly (d A) .poly (dT) 
and B^A in poly(dG).poly(dC). These results agree with 
the assumed similarity of the A^B transitions in ex- 
periment and in water drop simulations. It appears 
that, with hydration numbers close to stability of the 
B-form, these two fragments have qualitatively different 



counterion distributions, with a characteristic accumu- 
lation of Na"*" ions next to the opening of the minor 
groove in poly (dA). poly (dT). In fibers, these ion po- 
sitions can be additionally stabilized by intermolecular 
contacts, which would explain an increased persistence 
of poly (dA). poly (dT) DNA towards A-form in crystalline 
and amorphous fibers as compared to solution conditions. 



Results 

Comparative dynamics of A<^B transitions 

Fig. n] shows the sugar pucker dynamics in a series 
of representative trajectories of G12. Transitions from 
C2'-endo to C3'-endo are observed starting already from 
HN42. At first they occur only in the purine strand and 
for HN33 and HN29 an intermediate state is obtained 
with this strand almost entirely in the C3'-endo confor- 
mation. A complete transition to A-form is reached in 
the middle of the duplex starting from HN25. This tran- 
sition midway point is similar to that reported for the 
dodecamer CGCGAATTCGCG^. However, here the A- 
DNA backbone appears at much higher hydration num- 
bers and, as a result, the transition looks smoother and 
less cooperative. Moreover, the final A-DNA structure 
spreads from the center almost to entire fragment so 
that the B-philicity of the DNA ends, generally observed 
in experiments and calculation a^'^i^^ , is less significant. 
All this agrees well with the known A-philicity of the 
poly(dG).poly(dC) DNAS. 

Fig. ISlshows a titration-Hke pattern of these transitions 
as followed by different structural parameters. Different 
traces all exhibit S- or Z-shaped profiles sometimes with 
a very distinct transition zone^^. The A-DNA structures 
observed under low hydration are very close to the canon- 
ical conformation, with final RMSD values below 2 A. In 
contrast, under the highest hydration, an underwound 
B-DNA is observed with a strong bias towards A-form. 
Similar deviations were reported earher by others^SiffliSl 
and they are partially due to the known negative force- 
field bias of the Twisis^S,. The A^B transitions for G12 
are easily achieved in both directions and the points in 
Fig. 12 corresponding to the transition zone could be re- 
produced in alternative simulations starting from the A- 
DNA state. 

Similar results for A12 are shown in the next two fig- 
ures. The most striking difference from G12 could not 
be shown here, namely, that we failed to obtain a B— >A 
transition for this fragment in spite all efforts. At the 
same time, the A-form was perfectly stable under low 
hydration and, for HN21 and higher numbers, the A^B 
transitions passed without problems. Therefore, for Fig. 
El the trajectories chosen all started from the A-DNA 
conformation and we may add that, with the exception 
of HN13 and HN17, the B-DNA dynamics was also sta- 
ble and converged to the same structures. Comparison 
of Fig. El with Fig. readily reveals a stretch of hy- 
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FIG. 1: Dynamics of sugar pucker pseudorotatior^ in G12 under different hydration numbers. The boundary pucker values 
(ca 0°and 180°) are assigned the black and the white colors, respectively, with intermediate values mapped linearly to the 
gray scale levels. Each base pair step is characterized by a column consisting of two sub-columns, with the left sub-columns 
referring to the sequence written above in 5'-3' direction from left to right with the time averaged phases given on top. The 
right sub-columns refer to the complementary sequence shown below together with the corresponding time averaged phases. 
All trajectories shown in this figure started from the B-DNA state except the HN13. The last simulation started from the final 
state of HN17 because the original B-DNA state is unstable and prone to collapse under that low hydration. 



dration numbers were opposite transitions are observed 
in G12 and A12, that is the simulations reproduce the 
expected relative diflFerence between these two fragments 
as regards their preference towards A- vs B-DNA. 

Fig. exhibits a pattern qualitatively similar to that 
in Fig. 121 with a few exceptions. Again the computed 
A-DNA structures are very close to the canonical con- 
formation, but in this case the computed B-DNA are 
also relatively similar to the canonical one, even though 
the twist remains strongly underestimated. These traces 
were obtained by increasing hydration as shown in Fig. 
m In the opposite sense they could be also reproduced 
except for the hydration numbers below 21. In conditions 
HN17 and HN13, repetitive trajectories of A12 starting 
from B-DNA resulted in opening of terminal base pairs 



followed by irregular deformations that have been inter- 
preted as beginning of denaturation. Interestingly, when 
this limiting hydration is approached from the right in 
Fig. 01 the data show no trend towards the A-form. In- 
stead some parameters are nearly stable while the other 
even exhibit a trend in an opposite direction, suggest- 
ing that the B-DNA structure is progressively stabilized 
before the eventual base pair opening and collapse. 



Evolution of Na^ distributions around DNA 

The distributions of counterions around DNA are com- 
pared for G12 and A12 in Fig. The patterns shown 
here are interpreted as follows. The duplexes are entirely 
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FIG. 2: Quasi-static pattern of A^B transitions for G12 as 
monitored by different structural parameters. For each hy- 
dration number, a long MD trajectory was computed and the 
last 1 ns average structures were analyzed to obtain the nec- 
essary values. The two top plates show the average groove 
widths measured as described elsewhere—*—. Several helical 
parameters are all computed by Curves^, with the local val- 
ues for the Slide and global ones for the rest. The RMSD from 
the canonical A and B-DNA conformations is shown by cir- 
cles and squares, respectively. All distances are in angstroms 
and angles in degrees. The dotted lines indicate canonical A- 
and B-DNA levels. 



covered by water, therefore, the radial distribution of wa- 
ter oxygens shows the evolution of the available volume 
not occupied by DNA. Even though the relative sizes 
of water and Na"*" are not identical, these profiles give a 
rough estimate of free space that can be sampled by ions, 
and, in the absence of specific DNA-ion interactions, the 
solid and dotted traces in Fig. should have similar 



peak positions. In contrast, radial zones of strong Na"*"- 
DNA interaction produce strong separate peaks. Note 
also that the height of each peak in a cylindrical distribu- 
tion should be multiplied by the corresponding distance 
when their relative weights are estimated. 

The characteristic distributions for B-DNA are best 
seen in the several top plates of A12. It does not 
differ from similar distributions earlier reported for 
CGCGAATTCGCG^a. The three broad water peaks at 
approximately 4, 8 and 12.5 A correspond, respectively, 
to the first water layer in the major groove, the next few 
layers in both major and minor grooves and the bulk wa- 
ter outside DNA radius which is about 10 A for both A- 
and B-DNA. The first Na"*" peak at 5 A is a sum of sev- 
eral different contributions. It involves ions that interact 
with bases directly in depth of the minor groove as well 
as those in contact with N7 atoms of purine bases in the 
major groove. Mobile ions in the second hydration shell 
in the major groove also contribute to this peak because 
Na-OW distances are somewhat shorter than OW-OW 
ones. In cylindrical distributions in Fig. these qualita- 
tively different contributions could not be resolved and 
that is why Na"*" peaks at 5 A are found in B- as well 
as A-DNA distributions and often vary in height. The 
outer peak at 12.5 A includes highly mobile ions around 
DNA involved in non-specific phosphate screening. Its 
prominent part that becomes less mobile under low hy- 
dration corresponds to Na"^ positions along the minor 
groove between the two phosphate strands. Finally, an 
intermediate peak at approximately 8 A is generally pro- 
duced by mobile ions in the minor groove. Ion entering 
into the minor groove was relatively rare, perhaps be- 
cause the starting Na"*" distributions sampled only outer 
positions, which explains why such peaks are less well 
reproduced between the trajectories. 

The characteristic A-DNA distributions can be seen in 
the bottom plates of A12 as well as in many plates of 
G12. In this case the center of the DNA cylinder is filled 
with water and, consequently, accessible for ions. Most of 
them rest inside the DNA radius. As regards the relative 
weights, the most prominent peak is found at 10 A and 
it corresponds to Na"*" sandwiched between the opposite 
phosphate groups in the opening of the major groove. 
It is the accumulation of counterions at this region that 
causes the B— >A transition^fi. 

For G12 even under high hydration the Na"*" distribu- 
tion is rather different from typical B-DNA. These struc- 
tures have a strong negative Xdisp, therefore, the helical 
axis is shifted to the major groove and is accessible to 
water and ions. Even though in HN82 most of Na"'" stay 
outside the DNA radius their population inside the ma- 
jor groove is already significant. With hydration numbers 
reduced, the counterions are progressively pushed inside 
DNA. The Na"*" peak corresponding to the opening of the 
major groove emerges gradually and starting from HN25 
it becomes predominant. All this corresponds well to the 
relatively smooth B— >A transition characterized above 
in Fig. n] and Fig. [2 The A12 patterns in the right 
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FIG. 3: Dynamics of sugar pucker pseudorotation in A12 under different hydration numbers. All trajectories shown in this 
figure started from the A-DNA state . 



column are radically different. In this case both water 
and Na+ distributions retain the characteristic B-DNA 
shapes even with the hydration number reduced to 21. 
The ions are not pushed inside DNA when the outer wa- 
ter shell is reduced, and the Na"*" peak at 12.5 A remains 
dominant. 

An evident demonstration of the distinction between 
these two DNA fragments as regards their relative ca- 
pacities to adopt A- and B-forms is given in Fig. Elwhere 
the snapshots are shown for A12 and G12 at the end of 
dynamics simulated in drops of the same size (HN25). 
These structures and the counterion distributions con- 
firm the main aspects of dynamics illustrated by previous 
figures. For A12 we have a B-DNA with a narrow minor 
groove near both ends and a widening in the middle. A 



large number of the Na"^ ions stays outside DNA in front 
of the minor groove. In contrast, G12 gives a typical 
A-DNA conformation with a layer of Na"*" sandwiched 
between the opposite phosphate groups in the narrow 
major groove. This "electrostatic sandwich" provides the 
main driving force of the B— >A transition according to 
the condensation mechanisroiSiSfliSi. 



Persistence of poly(dA).poly(dT) DNA to A-form 

A series of additional calculations was carried out in or- 
der to check the origin of the A12 persistence with respect 
to A-DNA. It was possible that, the poly (dA). poly (dT) 
was intrinsically unstable with respect to base pair open- 
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FIG. 4: Quasi-static pattern of 
The notation is as in Fig. [2| 



>B transitions for A12. 



ing during the B^A transition due to mechanical stretch- 
ing of base pairs under intermediate hydration num- 
bers. In solution experiments, DNA commonly ex- 
hibits reduced thermostability under intermediate alco- 
hol concentrations^^. It can even denaturate and then 
renaturate when the water activity is further reducedS^. 
To probe this transitional instability several long trajec- 
tories were run for A-DNA close to the intermediate hy- 



FIG. 5: Cylindrical radial distribution functions for water 
oxygens (solid lines) and Na"*" ions (dashed lines) around G12 
and A12. DNA structures saved in the last nanosecond of 
dynamics together with surrounding water and counterions 
were superimposed with the canonical canonical B-DNA with 
the global coordinate OZ direction as the common helical axis. 
The Na"*" ions and water oxygens were counted in coaxial 0.1 
A thick cylinders. The distributions are volume normalized, 
that is scaled with a factor of 1/r, and the final plots were 
area normalized. The corresponding hydration numbers are 
shown on the right. 
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FIG. 6: Stereo snapshots of the final computed states of A12 and G12 obtained with the same drop size corresponding to the 
hydration number of 25. The Na"*" ions are shown by spheres. 



dration. Other simulations involved forced opening of 
one base pair to check if this could induce denaturation. 
In all these tests the A-DNA remained stable. 

The opening of terminal base pairs could also result 
from kinetic instability of B-DNA when it is suddenly 
transfered under low hydration. Such an opening is some- 
times observed even for G12 when trajectories start from 
B-DNA in HN13. To check this possibility, simulations 
were continued from the final B-DNA state of the HN21 
trajectory, with the amount of water reduced smoothly in 
several steps by removing each time only molecules not 
in contact with ions. This slow "evaporation" also led 
to opening of a terminal base pair under an intermediate 
hydration, with the structure remaining in the B-form. 

Finally, it was possible that the opening of terminal 



base pairs masks the B^A transition kineticahy, that is 
the molecule just cannot stay intact long enough for the 
transition to occur. To check this the two opposite ter- 
minal AT base pairs were replaced by GC and dynamics 
was simulated starting from B-DNA in HN17 and HN13 
water drops. In these simulations no base pair opening 
was observed. For HN17 two trajectories were computed 
of 11 and 7 ns, respectively, with different random initial 
distributions of counterions. In both cases the duplex 
remained firmly in the B-DNA state. For HN13 the tra- 
jectory was continued to 5.5 ns. During this time water 
started to leave the minor groove and the duplex suffered 
irregular deformations that drove it to a shrunk collapsed 
structure with partially closed grooves, and yet the tran- 
sition to the A-form did not occur. 
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The foregoing pattern is in a remarkable correspon- 
dence with the well documented poly (dA). poly (dT) pho- 
bicity towards the A-formS*2&. Unfortunately, this finding 
does not yet mean that a simple mechanism of this resis- 
tance can be readily extracted from simulations. It ap- 
peared that the B^A transitions in poly (dA). poly (dT) 
fragments become possible when the chain length is in- 
creased. In a 16-mer fragment of poly (dA). poly (dT) 
under HN17, opening of several terminal base pairs at 
both ends was also observed, however, this was accom- 
panied by a B^A transition in a few central residues. 
When the two terminal base pairs at the opposite ends 
were replaced with GC the resulting 16-mer fragment 
CGAAAAAAAAAAAAGC under the same hydration 
showed a normal B— >A transition in the middle 12-mer 
A-tract with the same general features as described be- 
fore. With the length of the A-tract further increased to 
18 base pairs, a B^A transition could also be obtained 
under HN21. 



Discussion 

The experimental studies of the B/A polymorphism 
in solution earlier revealed a surprisingly complex pat- 
tern of the double heHcal DNA behavior depending upon 
its sequence as well as environment conditions such as 
the water activity, the temperature, the type and the 
concentration of counterions. Perhaps not all effects ob- 
served in vitro are biologically relevant, but the excep- 
tional role of the DNA molecule makes any information 
concerning its properties potentially important. Unfor- 
tunately, a number of intriguing experimental observa- 
tions in this field could not be accurately interpreted in 
terms of DNA structure because the direct NMR and 
X-ray methods are not always applicable and often lim- 
ited in accuracy. Molecular dynamics simulations of the 
type described here can help to clarify some long stand- 
ing issues in this domain provided that the mechanism of 
A-<-»B transitions is similar in silico and in vitro. 

Under low hydration, dynamic DNA structure sys- 
tematically converges to the A-form conformation that 
is very similar for different sequences and very close to 
the canonical conformation observed in fiber crystals. In 
contrast, a much more significant and sequence depen- 
dent divergence from experimental structures is system- 
atically obtained for B-DNA. There is a surprising corre- 
spondence between this difference and experimental data. 
The high regularity and the absence of sequence effect for 
A-form was noticed long ago for X-ray fiber diffraction 
patterns^ and later confirmed in the ensemble of single 
crystal A-DNA structures^. Thus, in both experiment 
and simulations, the A-form of DNA is virtually insen- 
sitive to the base pair sequence. In simulations it ap- 
pears also to be relatively insensitive to the accuracy of 
the forcefield that do not allow MD trajectories to come 
closer to experimental B-DNA conformations. These ob- 
servations may have a common physical origin. Notably, 



they can be attributed to the A-form being dominated 
by the "electrostatic sandwich" in the major groove. 
The strong interactions of phosphate groups with metal 
cations effectively impose geometric constraints upon the 
inter-phosphate distances and suppress all other factors 
that might affect the overall structure. Similar types of 
interactions are also possible for the A-form stretches ob- 
served in complexes with proteins. In all such structures 
resolved by now the major groove of A-DNA is exposed 
to solvent whereas extensive protein DNA contacts are 
observed in the minor groove. 

The systematic deviations of the computed B-DNA 
conformations from the canonical one are similar to the 
recent simulations for analogous sequences^SiSSiSi. The 
Cornell et al parameters'^ have a well documented ten- 
dency to underestimate the average helical twisfeS^, and 
for the sequences studied here this bias is perhaps the 
largest. For poly(dG).poly(dC) this results in struc- 
tures with very A- like helical parameters. At the same 
time, sugar pseudorotation dynamics samples mainly 
from South/East phases characteristic of B-DNA (see 
Fig. nj. Earlier such structures were interpreted as 
strongly underwound B-DNA or as A-DNA with B-like 
backboneSSiSMi. Unfortunately, it is rather difficult to 
tell exactly how strong is the deviation of these conforma- 
tions from experimental data. They are rather different 
from G-tracts found in a large number of A-DNA single 
crystal structures. However, these structures are found in 
a unique A-DNA packing characterized by very special 
DNA-DNA interactions in wide minor grooves^. The 
same trend is not seen in short G-tracts found in protein- 
DNA complexes. For instance, one of the resolved struc- 
tures of the Y-family DNA polymerase contains a DNA 
fragment with an exposed GGGGG terminus'^. This G- 
tract features a normal B-DNA conformation with South 
sugar puckers and the average Twist of 34.2°rather close 
to that of generic B-DNA. Unfortunately, the last value 
is much larger than the measured solution value (ca 32.7° 
for long G-tracta^S,) suggesting that the DNA structure 
may again be affected by the crystal environment. On 
the other hand, solution NMR and Circular dichroism 
(CD) studies indicate that the G-tract duplexes com- 
puted with the Cornell et al force field resemble exper- 
imental structures^Siffi. The corresponding CD spectra 
are very similar to typical A-DNA suggesting that the 
high hydration conformation of poly(dG).poly(dC) al- 
ready has an A-like base pair stacking. At the same 
time, the sugar phases determined by NMR are all at 
the South. These conformations can hardly be assigned 
to the A-form because a cooperative transition to a gen- 
uine A-form is still distinguishable in solution as well as 
in crystalline fibers although the corresponding spectral 
changes are small^SiSMS. These subtle changes can well 
result from a transition pattern similar to that in Fig. 
IS with sugar phases switching to the North and grooves 
changing their width, but relatively minor shifts in helical 
parameters. 

The poly (dA). poly (dT) conformation at the high hy- 
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dration limit are much closer to the canonical B-DNA 
structure although, in agreement with earlier report^iii, 
it is also underwound by ca 4° with respect to the exper- 
imental solution value22i22iH. Our simulations managed 
to reproduce the well-documented relative A-form pho- 
bicity of this sequence in several different aspects. No- 
tably, there is a well defined range of hydration num- 
bers where poly (dA). poly (dT) and poly(dG).poly(dC) 
are stable in B- and A-forms, respectively. Dodecamer 
A12 resembles poly (dA). poly (dT) properties in crys- 
taUine fibers in that it refuses to go to the A-form and 
is prone to denaturation. At the same time, longer 
poly (dA). poly (dT) fragments with GC termini adopt the 
A-form without great problems, which agrees with exper- 
iments in solution--. It seems rather perplexing that a 
B— >A transition for A12 could not be obtained in spite of 
the fact that its A-form was found to be perfectly stable 
under low hydration. It is known, however, that in some 
experimental conditions, the B^A transition in DNA 
can go via temporary condensation or denaturation^. 
If the lowest energy transition pathway goes like that 
in our conditions, the molecular dynamics simulations 
would not be able to reach the final state because of the 
limited duration of trajectories. One should also keep in 
mind that the experimental B^A transitions in DNA re- 
portedly involve macroscopically slow relaxation phases 
and hysteresis eSectsS^^. Even though they usually con- 
cern long DNA chains, they may be due to local effects 
like slow ion diffusion especially under low hydration. 

Our results feature the complexity of the 
poly (dA). poly (dT) DNA as regards the B/A poly- 
morphism, but they leave unanswered the question 
concerning the main physical origin of such extraordi- 
nary behavior. A number of factors can contribute to 
the poly (dA). poly (dT) persistence with respect to the 
B— >A transition. The sequence dependent propensities 
towards A and B-forms should certainly originate from 
of the specific base stackingSiifliii. However, the variety 
and the character of the observed environmental effects 
suggest that subtle structural deviations prompted 
by specific base-base contacts are amplified by strong 
interactions with water and counterions. For instance, 
the narrow minor groove profile in the B-form of 
poly (d A) .poly (dT) may be sealed by the hydration spine 
structurei^. The same feature should produce a twofold 
effect upon the counterion distribution around DNA^iiii. 
The free positive ions should tend to accumulate in front 
of the narrow minor groove, and, simultaneously, their 
concentration in the major groove is reduced. Both these 
features should stabilize the B-form. Earlier experimen- 
tal and computational observations also suggest that a 
special role should be played by thymine methyl groups. 
They hinder a negative slide movement involved in the 
B— »A transitionii and form a continuous non-polar 
cluster in the major groove of B-DNA that should 
provide additional hydrophobic stabilization^^. The 
same feature effectively reduces the accessible volume of 
the major groove and can prevent accumulation of free 



solvent cationsSS. 

Many years ago, the contrasting propensities of 
poly (dA). poly (dT) and poly(dG).poly(dC) DNA to 
adopt A- and B-forms presented the first experimental 
demonstration of sequence dependent properties of the 
double helical DNA structure^i. Since then the repertory 
of reported sequence effects has many times increased, 
and yet the exact physical origin of this particular differ- 
ence remains controversial. It is shown here that A^B 
transitions observed in water drop simulations exhibit 
clear trends qualitatively similar to the long known ex- 
perimental observations. These results corroborate the 
putative general role of the intra-duplex electrostatic con- 
densation mechanism for A-f->B transitions in DNA in 
vitro and suggest that future studies in the same direc- 
tion can give more definite answers to the issues discussed 
here. 



Methods 

A series of MD simulations was carried out for dou- 
ble helical dodecamer fragments poly (dA). poly (dT) and 
poly(dG).poly(dC) referred to as A12 and G12, respec- 
tively, as well as a few derivatives of these sequences. 
Water drops of 300, 400, 500, 600, 700, 800, 1000, and 
2000 molecules were used, which gives approximate hy- 
dration numbers (water/nucleotide) 13, 17, 21, 25, 29, 
33, 42, and 83, respectively. Additional simulations with 
a hydration number around 171 were carried out by us- 
ing periodical boundary conditions. Below for brevity 
the simulation conditions outlined above are referred to 
as HN13, HN17, and so fourth. All calculations were 
continued long enough to obtain convergent average pa- 
rameters of the final states. In this way a quasi-static 
pattern of B<-^A transitions can be reproduced resem- 
bling in vitro titration experiments 

The simulation protocols were similar to the earlier 
water drop simulations^SiS^. We use the internal coordi- 
nate molecular dynamics (ICMD) methodS*2^ adapted 
for DN ATTiTff with the time step of 0.01 ps. In this ap- 
proach, the DNA molecule has all bond length and al- 
most all bond angles fixed at their standard values. The 
only variable bond angles are those centered at the sugar 
C1',.--,C4', and 04' atoms, which assures the fiexibil- 
ity of the furanose rings. In contrast, bases, thymine 
methyls, and phosphate groups move as articulated rigid 
bodies, with only rotations around single bonds allowed. 
The highest frequencies in thus obtained models are ad- 
ditionally balanced by increasing rotational inertia of 
the lightest rigid bodies as described earlies22i£&. The 
possible physical effects of the above modifications have 
been discussed elsewher oT^ii^^ . The electrostatic inter- 
actions are treated with the recent version of the SPME 
methodSi specifically adapted for infinite vacuum bound- 
ary conditions'^. The common values of Ewald parame- 
ters were used, that is 9 A truncation for the real space 
sum and /3 ~ 0.35. Reference simulations with periodical 
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boundaries were carried out as described beforeS^, with 
the standard SPME method in NVT ensemble conditions 
with a rectangular unit cell of 45 x 45 x 65 A under normal 
water density. 

The initial states were prepared as follows. The canon- 
ical A and B-DNA^i were used as standard initial con- 
formations. The DNA molecule was first immersed in a 
large rectangular TIP3B^ water box of and next exter- 
nal solvent molecules were removed by using a spherical 
distance cut-off from DNA atoms. The cut-off radius was 
adjusted to obtain the desired number of water molecules 
remaining. The drop was neutralized by randomly plac- 
ing Na"*" ions at water positions selected so that their 
distances from DNA were 5 A or larger. The initial coun- 
terion distribution was pre-equilibrated by running 1 ns 
dynamics in water drops of 500 molecules for A-DNA and 
800 molecules for B-DNA, with DNA atoms weakly re- 
strained to their initial positions. The final drop size was 
adjusted by adding or removing water from the surface. 



The described procedure was intended to assure the start 
of dynamics from closely similar states regardless of the 
drop size. 

Every system was energy minimized first with the so- 
lute held rigid and then with all degrees of freedom. 
Dynamics were initiated with the Maxwell distribution 
of generalized momenta at low temperature. The sys- 
tem was next slowly heated to 250 K and equilibrated 
during several picoseconds. Production trajectories were 
computed with the temperature bound to 293 K by the 
Berendsen algorithn:>24 with a relaxation time of 10 ps. 
For better comparison with earlier simulations of A<->B 
transitions, the original Cornell et al. force filedS^ was 
used. Duration of production runs varied from 2 to 25 
ns depending upon the observed character of dynamics. 
The conformations were saved with a 2.5 ps interval. 
The results were analyzed with in-house routines and the 
Curves program^. 



^ Saenger, W. (1984). Principles of Nucleic Acid Structure, 

Springer- Verlag, New York. 
^ Ivanov, V. I. & Minchenkova, L. E. (1995). The A-form of 

DNA: in search of biological role (a review). Mol. Biol. 28, 

780-788. 

^ Lu, X. J., Shakked, Z., & Olson, W. K. (2000). A-form 
conformational motifs in ligand-bound DNA structures. J. 
Mol. Biol. 300, 819-840. 

^ Franklin, R. E. & Gosling, R. G. (1953). Molecular config- 
uration in sodium thymonucleate. Nature 171, 740-741. 

^ Leslie, A. G. W., Arnott, S., Chandrasekaran, R., & 
Ratliff, R. L. (1980). Polymorphism of DNA double he- 
lices. J. Mol. Biol. 143, 49-72. 

Piskur, J. & Rupprecht, A. (1995). Aggregated DNA in 
ethanol solution. FEBS Lett. 375, 174-178. 
Tunis-Schneider, M. J. & Maestre, M. F. (1970). Circular 
dichroism spectra of oriented and unoriented deoxyribonu- 
cleic acid films: A preliminary study. J. Mol. Biol. 52, 
521-541. 

* Ivanov, V. I., Minchenkova, L. E., Minyat, E. E., Frank- 
Kametetskii, M. D., & Schyolkina, A. K. (1974). The B to 
A transition of DNA in solution. J. Mol. Biol. 87, 817-833. 

^ Calladine, C. R. & Drew, H. R. (1984). A base-centered 
explanation of the B-to-A transition in DNA. J. Mol. Biol. 
178, 773-782. 

Mazur, J., Sarai, A., & Jernigan, R. L. (1989). Sequence 
dependence of the B-A conformational transition of DNA. 
Biopolymers 28, 1223-1233. 

Hunter, C. A. (1993). Sequence-dependent DNA structure. 
The role of base stacking interactions. J. Mol. Biol. 230, 
1025-1054. 

Ivanov, V. I., Zhurkin, V. B., Zavriev, S. K., Lysov, Y. P., 
Minchenkova, L. E., Minyat, E. E., Frank-Kamenetskii, 
M. D., & Schyolkina, A. K. (1979). Conformational possi- 
bilities of double-helical nucleic acids: Theory and experi- 
ment. Int. J. Quantum Chem. 16, 189-201. 
" Cheatham, III, T. E. & Kollman, P. A. (1996). Observation 
of the A-DNA to B-DNA transition during unrestrained 
molecular dynamics in aqueous solution. J. Mol. Biol. 259, 



434-444. 

" Yang, L. & Pettitt, B. M. (1996). B to A transition of 
DNA on the nanosecond time scale. J. Phys. Chem. B 100, 
2564-2566. 

Cheatham, III, T. E., Crowley, M. F., Fox, T., & Kollman, 
P. A. (1997). A molecular level picture of the stabilization 
of A-DNA in mixed ethanol-water solutions. Proc. Natl. 
Acad. Sci. U.S.A. 94, 9626-9630. 
1" Cheatham, III, T. E. & Kollman, P. A. (1997). Insight 
into the stabilization of A-DNA by specific ion association: 
Spontaneous B-DNA to A-DNA transitions observed in 
molecular dynamics simulations of d[ACCCGCGGGT]2 in 
the presence of hexaamminecobalt(III). Structure 5, 1297- 
1311. 

" Cieplak, P., Cheatham, III, T. E., & Kollman, P. A. 

(1997) . Molecular dynamics simulations find that 3'- 
phosphoramidate modified DNA duplexes undergo B to A 
transition and normal DNA duplexes an A to B transition. 
J. Am. Chem. Soc. 119(29), 6722-6730. 

Jayaram, B., Sprous, D., Young, M. A., & Beveridge, D. L. 

(1998) . Free energy analysis of the conformational prefer- 
ences of A and B forms of DNA in solution. J. Am. Chem. 
Soc. 120, 10629-10633. 

" Sprous, D., Young, M. A., & Beveridge, D. L. (1998). 
Molecular dynamics studies of the conformational prefer- 
ences of a DNA double helix in water and an ethanol/water 
mixture: Theoretical considerations of the A-B transition. 
J. Phys. Chem. B 102, 4658-4667. 

Mazur, A. K. (2003). Titration in silico of reversible B^A 
transitions in DNA. J. Am. Chem. Soc. 125(26), 7849- 
7859. 

Feig, M. & Pettitt, B. M. (1998). Structural equilibrium 
of DNA represented with different force fields. Biophys. J. 
75, 134-149. 

Mazur, A. K. (2002). DNA dynamics in a water drop 
without counterions. J. Am. Chem. Soc. 124(49), 14707- 
14715. 

Feuerstein, B. G., Pattabiraman, N., & Marton, L. J. 
(1986). Spermine-DNA interactions: a theoretical study. 



11 



Proc. Natl. Acad. Sci. U. S. A. 83, 5948-52. 

Rouzina, I. & Bloomfield, V. A. (1998). DNA bending by 

small, mobile multivalent cations. Biophys. J. 74, 3152- 

3164. 

Xu, Q., Shoemaker, R. K., & Braunlin, W. H. (1993). In- 
duction of B- A transitions of deoxyoligonucleotides by mul- 
tivalent cations in dilute aqueous solution. Biophys. J. 65, 
1039-1049. 

2" Saenger, W., Hunter, W. N., & Kennard, O. (1986). DNA 
conformation is determined by economics in the hydration 
of phosphate groups. Nature 324, 385-388. 
Ivanov, V. I., Minchenkova, L. E., Schyolkina, A. K., & 
Poletaev, A. I. (1973). Different conformations of double- 
stranded nucleic acids in solution as revealed by circular 
dichroism. Biopolymers 12, 89-110. 

Zimmerman, S. B. (1983). The three-dimensional structure 
of DNA. Ann. Rev. Biochem. 51, 395-427. 
Becker, M. M. & Wang, Z. (1989). B->A transitions within 
a 5 S ribosomal RNA gene are highly sequence-specific. J. 
Biol. Ghem. 265, 4163-4167. 
^° Basham, B., Schroth, G. P., & Ho, P. S. (1995). An A- 
DNA triplet code: Thermodynamic rules for predicting A- 
and B-DNA. Proc. Natl. Acad. Sci. U.S.A. 92, 6464-6468. 
Tolstorukov, M. Y., Ivanov, V. I., Malenkov, G. G., Jerni- 
gan, R. L., & Zhurkin, V. B. (2001). Sequence-dependent 
B-A transition in DNA evaluated with dimeric and trimeric 
scales. Biophys. J. 81, 3409-3421. 

Peck, L. J. & Wang, J. C. (1981). Sequence dependence of 
the helical repeat of DNA in solution. Nature 292, 375- 
378. 

Rhodes, D. & Klug, A. (1981). Sequence-dependent helical 
periodicity of DNA. Nature 292, 378-380. 
Strauss, P., Gaillard, C, & Prunell, A. (1981). Heli- 
cal periodicity of DNA, poly (dA). poly (dT) and poly(dA- 
dT).poly(dA-dT) in solution. Eur. J. Biochem 118, 215- 
222. 

Hogan, M., LeGrange, J., & Austin, B. (1983). Dependence 
of DNA helix fiexibility on base composition. Nature 304, 
752-754. 

Nishimura, Y., Torigoe, C., & Tsuboi, M. (1986). Salt in- 
duced B to A transition of poly(dG)-poly(dC) and the sta- 
bilization of A form by its methylation. Nucleic Acids Res. 
14, 2737-2749. 

^'^ Behling, R. W. & Kearns, D. R. (1986). IH two- 
dimensional nuclear overhauser effect and relaxation stud- 
ies of poly (dA). poly (dT). Biochemistry 25, 3335-3346. 
Benevides, J. M., Wang, A. H., Rich, A., Kyogoku, Y., 
van der Marel, G. A., van Boom, J. H., & Thomas, 
Jr., G. J. (1986). Raman spectra of single crystals of 
r(GCG)d(CGC) and d(CCCGGG) as models for A DNA, 
their structure transitions in aqueous solution, and com- 
parison with double-helical poly (dG). poly (dC). Biochem- 
istry 25, 41-50. 

Sarma, M. H., Gupta, G., & Sarma, R. H. (1986). 500-MHz 
IH NMR study of poly (dG). poly (dC) in solution using one- 
dimensional nuclear Overhauser effect. Biochemistry 25, 
3659-3665. 

'*° Peticolas, W. L., Wang, Y., & Thomas, G. A. (1988). 
Some rules for predicting the base-sequence dependence 
of DNA conformation. Proc. Natl. Acad. Sci. U.S.A. 85, 
2579-2583. 

Buckin, V. A., Kankiya, B. I., Bulichov, N. V., Lebedev, 
A. v., Gukovsky, I. Y., Chuprina, V. P., Sarvazyan, A. P., 
& Williams, A. R. (1989). Measurement of anomalously 



high hydration of (dA)n.(dT)n double helices in dilute so- 
lution. Nature 340, 321-322. 

Buckin, V. A., Kankiya, B. I., Sarvazyan, A. P., & Uedaira, 

H. (1989). Acoustical investigation of poly (dA). poly (dT), 
[poly[d(A-T)], poly (A). poly (U) and DNA hydration in di- 
lute aqueous solutions. Nucleic Acids Res. 17, 4189-4203. 

*^ Brahms, S., Pritsch, V., Brahms, J. G., & Westhof, 
E. (1992). Investigations on the dynamics structures of 
adenine- and thymine-containing DNA. J. Mol. Biol. 223, 
455-476. 

Chalikian, T. V., Plum, G. E., Sarvazyan, A. P., 
& Breslauer, K. J. (1994). Infiuence of drug bind- 
ing on DNA hydration: acoustic and densimet- 
ric characterizations of netropsin binding to the 
poly (dAdT). poly (dAdT) and poly (dA). poly (dT) du- 
plexes and poly (dT). poly (dA). poly (dT) triplex at 25 
degrees C. Biochemistry 26, 8629-8640. 
Vorlickova, M., Subriana, J. A., Chladkova, J., Tejralova, 

I. , Huynh-Dinh, T., Arnold, L., & Kypr, J. (1986). 
Comparison of the solution and crystal conformations of 
(G+C)-rich fragments of DNA. Biophys J. 71, 1530-1538. 

*^ Ivanov, V. I., Minchenkova, L. E., Burckhardt, G., Birch- 
Hirschfeld, E., Pritzsche, H., & Zimmer, C. (1996). The 
detection of B-form/A-form junction in a deoxyribonu- 
cleotide duplex. Biophys. J. 71, 3344-3349. 
Hud, N. V. & Plavec, J. (2003). A unified model for the 
origin of DNA sequence-directed curvature. Biopolymers 
69, 144-159. 

Hud, N. V. & Polak, M. (2001). DNA-cation interactions: 
The major and minor grooves are flexible ionophores. Curr. 
Opin. Struct. Biol. 11, 293-301. 

*^ Dechering, K. J., Cuelenaere, K., Konings, R. N., & Le- 
unissen, J. A. (1998). Distinct frequence-distributions of 
homopolymeric DNA tracts in different genomes. Nucleic 
Acids Res. 26, 4056-4062. 

'° Trifonov, E. N. & Sussman, J. L. (1980). The pitch of chro- 
matin DNA is reflected in its nucleotide sequence. Proc. 
Natl. Acad. Sci. U.S.A. 77(7), 3816-3820. 

®^ Widlund, H. R., Cao, H., Simonsson, S., Magnusson, E., 
Simonsson, T., Nielsen, P. E., Kahn, J. D., Crothers, D. D., 
& Kubista, M. (1997). Identification and characterization 
of genomic nucleosome-positioning sequences. J. Mol. Biol. 
267, 807-817. 

Vashakidze, R. P. & Prangishvili, D. A. (1987). Sim- 
ple repetitive sequences in the genomes of archaebacteria. 
FEBS Lett. 216, 217-220. 
®^ Altona, C. & Sundaralingam, M. (1972). Conformational 
analysis of the sugar ring in nucleosides and nucleotides. 
A new description using the concept of pseudorotation. J. 
Am. Chem. Soc. 94(23), 8205-8212. 

Mazur, A. K. (1999). Internal correlations in minor groove 
profiles of experimental and computed B-DNA conforma- 
tions. J. Mol. Biol. 290, 373-377. 

Lavery, R. & Sklenar, H. (1988). The definition of general- 
ized helicoidal parameters and of axis curvature for irreg- 
ular nucleic acids. J. Biomol. Struct. Dyn. 6(1), 63-91. 
^'^ Minchenkova, L. E., Schyolkina, A. K., Chernov, B. K., & 
Ivanov, V. I. (1986). CC/GG contacts facilitate the B to 
A transition of DNA in solution. J. Biomol. Struct. Dyn. 
4, 463-476. 

The current AMBER parameters reportedly overestimate 
the stability of cytosine C2'-endo sugar pucker with re- 
spect to the C3'-endo^. This forcefield defect is perhaps 
responsible for a delayed transition of cytosine sugars to 



12 



C3'-endo and it may somewhat reduce the A-philicity of 
poly(dG).poly(dC) DNA. 

The quasi-static profiles in our calculations should be dis- 
tinguished from familiar S-shaped experimental plots of co- 
operative A^B transitionaSi. In the latter case they result 
from two different contributions: (i) the shift in relative 
populations of A- and B-forms in the ensemble of DNA 
conformations, and (ii) small deformations of the A and 
B-conformations under varied hydration. Our simulations 
probe only the second contribution. Its relative weight in 
experiments is not known well. The Ising model of coopera- 
tive A^B transitions considers only the first contributiortS-. 
The second one is effectively neglected by assuming that in 
small DNA fragments the transition is described by a step 
function. It is known, however, that, beyond the transition 
zone of hydration numbers, conformations of both B- and 
A-DNA change in a smooth and non-cooperative manner 
and that these changes are rather significant^^J—. Evalua- 
tion of equilibrium populations of A and B-forms from this 
types of simulations rests beyond the current possibilities, 
nevertheless, reversion of A<->B transitions in a single MD 
trajectory can be observed in very long simulations under 
intermediate hydration (unpublished results). 
Trantirek, L., Stefi, R., Vorlickova, M., Koca, J., Sklenar, 
v., & Kypr, J. (2000). An A-type double helix of DNA 
having B-type puckering of the deoxyribose rings. J. Mol. 
Biol. 297, 907-22. 

Stefi, R., Trantirek, L., Vorlickova, M., Koca, J., Sklenar, 
v., & Kypr, J. (2001). A-like guanine-guanine stacking in 
the aqueous DNA duplex of d(GGGGCCCC). J. Mol. Biol. 
307, 513-24. 

Lankas, F., Cheatham, III, T. E., Spaskova, N., Hobza, P., 
Langowskii, J., & Sponer, J. (2002). Critical effect of the 
N2 amino group on structure, dynamics, and elasticity of 
DNA plypurine tracts. Biophys. J. 82, 2592-2609. 
Cheatham, III, T. E., Cieplak, P., & KoUman, P. A. (1999). 
A modified version of the Cornell et al. force field with im- 
proved sugar pucker phases and helical repeat. J. Biomol. 
Struct. Dyn. 16, 845-862. 

Ivanov, V. I., Krylov, D. Y., & Minyat, E. E. (1985). Three 
state diagram for DNA. J. Biomol. Struct. Dyn. 3, 43-55. 
Vorlickova, M., Sagi, J., Hejtmankova, I., & Kypr, J. 
(1991). Alkyl substituent in place of the thymine methyl 
group controls the A-X conformational bimorphism in 
poly(dA-dT). J. Biomol. Struct. Dyn. 9, 571-8. 
Suzuki, M., Amano, N., Kakinuma, J., & Tateno, M. 
(1997). Use of a 3D structure data base for understand- 
ing sequence-dependent conformational aspects of DNA. 
J. Mol. Biol. 274, 421-35. 

Cornell, W. D., Cieplak, P., Bayly, C. I., Gould, I. R., 
Merz, K. M., Ferguson, D. M., Spellmeyer, D. C, Fox, 
T., Caldwell, J. W., & Kollman, P. A. (1995). A second 
generation force field for the simulation of proteins, nucleic 
acids and organic molecules. J. Am. Chem. Soc. 117(19), 
5179-5197. 

Wahl, M. C. & Sundaralingam, M. (1999). A-DNA du- 
plexes in the crystal. In Stephen Neidle, (ed.), Oxford 
Handbook of Nucleic Acid Structure, pp. 117-144 Oxford 
University Press New York. 

Ling, H., Boudsocq, F., Woodgate, R., & Yang, W. (2001). 
Crystal structure of a Y-family DNA polymerase in action: 
a mechanism for error-prone and lesion-bypass replication. 



Cell. 107, 91-102. 

Arnott, S. & Seising, E. (1974). Letter: The structure of 
polydeoxyguanylic acid with polydeoxycytidylic acid. J. 
Mol. Biol. 88, 551-2. 

McConnell, K. J. & Beveridge, D. L. (2001). Molecular 
dynamics simulations of B'-DNA: Sequence effects on A- 
tract-induced bending and fiexibility. J. Mol. Biol. 314, 
23-40. 

Falk, M., Hartman, Jr., K. A., & Lord, R. C. (1962). Hy- 
dration of deoxyribonucleic acid. I. A gravimetric study. J. 
Am. Chem. Soc. 84, 3843-3846. 

Lindsay, S. M., Lee, S. A., Powell, J. W., Weidlich, T., 
Demarco, C, Lewen, G. D., Tao, N. J., & Rupprecht, A. 
(1988). The origin of the A to B transition in DNA fibers 
and films. Biopolymers 27, 1015-1043. 
Dickerson, R. E., Drew, H. R., Conner, B. N., Wing, R. M., 
Fratini, A. V., & Kopka, M. L. (1982). The anatomy of A-, 
B-, and Z-DNA. Science 216, 475-485. 
Pilet, J., Blicharski, J., & Brahms, J. (1975). Conforma- 
tions and structural transitions in polydeoxynucleotides. 
Biochemistry 14, 1869-76. 

Mazur, A. K. (1997). Quasi-Hamiltonian equations of mo- 
tion for internal coordinate molecular dynamics of poly- 
mers. J. Comput. Chem. 18(11), 1354-1364. 
Mazur, A. K. (2001). Internal coordinate simulation 
method. In Oren M. Becker, Alexander D. MacKerell, Jr., 
Benoit Roux, & Masakatsu Watanabe, (ed.). Computa- 
tional Biochemistry and Biophysics, pp. 115-131 Marcel 
Dekker New York. 

Mazur, A. K. (1998). Accurate DNA dynamics without 
accurate long range electrostatics. J. Am. Chem. Soc. 
120(42), 10928-10937. 

Mazur, A. K. (1999). Symplectic integration of closed chain 
rigid body dynamics with internal coordinate equations of 
motion. J. Chem. Phys. 111(4), 1407-1414. 
Mazur, A. K. (1998). Hierarchy of fast motions in protein 
dynamics. J. Phys. Chem. B 102(2), 473-479. 
Mazur, A. K., Sumpter, B. G., & Noid, D. W. (2001). 
Internal coordinate phase space analysis of macromolecular 
systems. Comput. Theor. Polym. Sci. 11(1), 35-47. 
Essmann, U., Perera, L., Berkowitz, M. L., Darden, T., 
Lee, H., & Pedersen, L. G. (1995). A smooth particle mesh 
Ewald method. J. Chem. Phys. 103, 8577-8593. 
Arnott, S. & Hukins, D. W. L. (1972). Optimised parame- 
ters for A-DNA and B-DNA. Biochem. Biophys. Res. Com- 
mun. 47, 1504-1509. 

Jorgensen, W. L., Chandreskhar, J., Madura, J. D., Im- 
pey, R. W., &i Klein, M. L. (1983). Comparison of simple 
potential functions for simulating liquid water. J. Chem. 
Phys 79, 926-935. 

Berendsen, H. J. C, Postma, J. P. M., van Gunsteren, 
W. F., DiNola, A., & Haak, J. R. (1984). Molecular dy- 
namics with coupling to an external bath. J. Chem. Phys. 
81, 3684-3690. 

Lee, C.-H., Mizusawa, H., & Kakefuda, T. (1981). Unwind- 
ing of double-stranded DNA helix by dehydration. Proc. 
Natl. Acad. Sci. U.S.A. 78, 2838-2842. 
Harmouchi, M., Albiser, G., & Premilat, S. (1990). 
Changes of hydration during conformational transitions of 
DNA. Eur. Biophys. J. 19, 87-92. 



